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Integrated computational materials engineering (ICME) is a useful approach for tailoring 
the performance of a material. For fiber-reinforced composites, not only do the properties 
of the constituents of the composite affect the performance, but so does the architecture 
(or microstructure) of the constituents. The generalized method of cells is demonstrated 
to be a viable micromechanics tool for determining the effects of the microstructure on 
the performance of laminates. The micromechanics is used to predict the inputs for a 
macroscale model for a variety of different fiber volume fractions, and fiber architectures. 
Using this technique, the material performance can be tailored for specific applications by 
judicious selection of constituents, volume fraction, and architectural arrangement given a 
particular manufacturing scenario. 
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MMCDM damage stiffness parameters 

MMCDM compressive damage stiffness parameters 

MMCDM tensile damage stiffness parameters 

Strain concentration matrix 

Subcell strain concentration submatrix 

Global RUC stiffness, Pa 

Normal compressive damage weighting parameters 
Tensie compressive damage weighting parameters 
Damage weighting parameters 
Subcell stiffness, Pa 
RUC dimensions, m 

MMCDM compressive scalar damage variables 

MMCDM tensile scalar damage variables 

MMCDM scalar damage variables 

Subcell dimensions, m 

MMCDM degradation factors 

MMCDM damage increment 

Damage strain increment 

Young’s modulus of matrix, Pa 

Axial and transvese moduli of fiber, Pa 

Undamaged Young’s moduli, Pa 
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Young’s moduli, Pa 

Undamaged transverse modulus of lamina, Pa 
Transverse modulus of lamina, Pa 

Tensile, compressive, and shear ST microdamage functions 

Thermodynamic driving force 

Shear modulus of fiber, Pa 

Undamaged shear moduli, Pa 

Shear moduli, Pa 

Instantaneous tangent stiffness, Pa 
Normalized tangent stiffness 
Total number of subcells 
Shear strain allowables 

ST microdamage ISV/dissipated potential, Pa 

Rate ST microdamage ISV/dissipated potential, Pa/s 

ST reduced microdamage ISV, Pas 

Elastic strain energy density, Pa 

Total work potential, Pa 

Compressive normal strain allowables 

Tensile normal strain allowables 

Normal strain allowables 

Macroscale composite coordinates 

Microscale RUC coordinates 

Subcell indices 

Engineering shear strain components 

Average global strains 

Vector containing all Subcell strains 

Subcell strains 

Damage strains 

Normal strain components 

Poisson’s ratio of matrix 

Axial and transverse Poisson’s ratios of fiber 

Undamaged Poisson’s ratios 

Poisson’s ratios 

Damage stress, Pa 

Subcell stresses Pa 

Normal stress components, Pa 


I. Introduction 

With the increased emphasis on reducing the cost and time to market of new materials, Integrated 
Computational Materials Engineering (ICME) has become a fast growing discipline within materials science 
and engineering. ICME is an integrated approach to the design of products and the materials which comprise 
them by linking material models at multiple time and length scales; such that manufacturing processes, which 
produce internal material structures that in turn influence material properties and allowables, can be tailored 
(engineered) to specific industrial applications. Figure 1 illustrates the interconnection of these scales and 
their cause/effect relationships; e.g., processing conditions produce a particular microstructure from which 
properties are obtained, which then dictate a specific structural performance. Note that the evolution of 
elliptical line types (i.e., dotted to dashed to solid line) are purposely included to visually suggest the 
level of maturity /understanding (from immature, to semi-mature, to mature, respectively) of modeling at 
each level of scale (both temporal and geometric). Furthermore, the figure demonstrates the importance 
of understanding the input and output at each scale in order to determine meaningful properties that are 
ultimately required by a structural analyst. Equally important is the fact that experiments (whether virtual 
or laboratory) performed at a given level can be viewed either as exploration or characterization experiments 
used to obtain the necessary model features or parameters, respectively, operating at the present and/or next 
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higher level, or, they can be used to validate the modeling methods employed to transition from the lower 
level to the given level. 



L«El£tb 

Figure 1: Schematic illustrating the multiple spatial and temporal scales active within material and structure. 


While there is a clear indication that ICME is growing, the realization of the successful implementation 
of ICME in daily work of researchers and engineers in industry and academia is still lacking. In this paper we 
will demonstrate how the general, synergistic, multiscale-modeling framework for composites (developed by 
the NASA Glenn Research Center (GRC) and known as MAC/GMC and FEAMAC, see 1 ) can be effectively 
utilized to link the material microstructure (e.g., constituent phase properties, volume fraction, fiber packing, 
etc.) to ply/laminate properties (mesoscale) and finally to performance (at the macroscale), see Figure 2, in 
an efficient and accurate manner to enable ”fit-for-purpose” tailoring of the composite material. The ability 
to localize and homogenize between scales with efficiency and accuracy make MAC/GMC and FEAMAC 
ideal candidates for ICME simulations in a multiscale environment in which the microstructure can be 
optimized spatially based on the local loading and environmental history. 



Figure 2: Illustration of relevant levels of scales for multiscale composite analysis. 


One promising lamina- level theory for modeling matrix microdamage is Schapery theory (ST), 2 ' 7 a 
thermodynamically-based work potential model. ST was successfully used previously to accurately predict 
the behavior of numerous off-axis unidirectional, angle-ply and multi-angle laminates. 3 4 However, lacking 
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explicit microstructural details such as fiber volume fraction and fiber packing, this theory cannot be utlized 
to design the material (i.e., for ICME) since any such modification would constitute a ’’new material” and 
require re-characterization of the model. Therefore, the scope of its usage is limited to strictly analysis. 

Micromechanics differs from macromechanics in that it enables one to account explicitly for variations 
in constituent material properties enhancements as well as microstructural architectural effects such as fiber 
volume content, fiber packing, void/defect content and orientations making micromechanics based models 
particularly attractive for use in ICME of composite materials. In this work, MAC/GMC, 8,9 which uses the 
generalized method of cells (GMC) micromechanics theory, is presented as a computational tool suitable for 
achieving the required microstructure to property and property to performance relationships necessary to 
enable the application specific tailoring of the carbon fiber-reinforced, polymer matrix composites (PMCs). 

Due to prior success in using ST, the current objective is to demonstrate that MAC/GMC can give 
comparable accuracy in predicting composite behavior while still enabling microstructural tailoring of the 
material. Since ST was formulated assuming, transversely isotropic lamina, it is not suitable for inclusion 
in a micromechanics model where the matrix phase requires a full 3-D, and often isotropic, formulation. 
Instead, a multi-axial mixed-mode continuum damage model 10 (MMCDM) is utilized to model the non- 
linear behavior of the matrix, due to the accumulation of microdamage within a PMC repeating unit cell 
(RUC). This micromechanics model is calibrated to match the same experimental data, which was used to 
characterize the microdamage input functions for ST. Results from the calibrated micromechanics model, 
using MMCDM, are compared to the ST and experimental results for additional off-axis, angle ply, and 
multi- angle laminates. 

Finally, manufacturing sensitive parameters such as fiber volume fraction and fiber packing are varied to 
observe the influence of these parameters on the ST input (microdamage function) and the performance of 
the multi-angle laminates. The analyses performed herein only considers matrix microdamage are invalid 
once meso-cracking (transverse cracking), fiber breakage, delamination or any other more severe damage 
mechanisms arise. 


II. The Generalized Method of Cells 


GMC, first developed in Ref. 11, offers a very computational efficient means of calculating the local 
fields within a composite RUC. It is completely general, and can admit any number of different phases with 
arbitrary deformation and damage constitutive laws, making it ideal for ICME applications. GMC has been 
utilized thoroughly for progressive damage and failure analysis (PDFA) of composite materials. 10, 12,13 As 
such, GMC is incredibly amenable for multiscale analysis of the same sort. 7, 14 16 A brief summary of the 
GMC fundamentals are presented in this section; the reader is referred to Ref. 1 for more details on GMC 
and the high-fidelity generalized method of cells (HFGMC). 

With GMC, an RUC must be identified within the composite microstructure, and the RUC is discretized 
into a number of subcells, see Figure 3. Each subcell contains a single constituent (obeying the constitutive 
law of that constituent), but there is no limit on the number of different constituents that may be present 
in the RUC. Displacement and traction continuity is enforced in an average, or integral sense at each of the 
subcell interfaces and the periodic boundaries of the RUC. These continuity conditions are used to formulate 
a strain concentration matrix A, 

e s = Ae (1) 

which gives all the local subcell strains (es) in terms of the six, global, average, applied strains e. 

Once the strain concentration matrix is formulated and the local subcell strains have been calculated, 
the local subcell stresses (cr^/U)) can b e calculated using the local constitutive law and the local subcell 
strains. 

Q.iaP'y) _ cO/U^O/U) (2) 

where <a, /?, and 7 are the subcell indices in the local 7 / 1 - , 7 / 2 -, and 7 / 3 -direct ions, respectively, C^ af3 ^ is the 
local subcell stiffness and contains the six local subcell strains. Eqs. (1) and (2) can be modified to 

admit inelasticity, however this manuscript focuses on stiffness reduction damage only. Finally, the global 
RUC stiffness can be calculated. 
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(a) Triply-periodic composite mi- (b) Method of cells idealization of 
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Figure 3 : Composite with repeating microstructure and arbitrary constituents. 


where N a , Np, and TV 7 are the total number of subcells in the yi -, 7/2-, and t/ 3 -directions, respectively, D , 
iJ, and L are the RUC dimensions in the 7/1 -, t/ 2 -, and t/ 3 -direct ions, d a , hp , and / 7 are the local subcell 
dimensions in the 7/1-, t/ 2 -, and 7/3-directions for subcell a/?7, and A^ 7 ) is a submatrix of the strain 
concentration matrix A corresponding to subcell 

III. Multiaxial Mixed-mode Continuum Damage Model for Matrix 

Constituents 


III. A. Damage Initiation 

The MMCDM, developed in Ref. 10 was implemented within MAC/GMC, and used to model progressive 
damage in the matrix constituent subcells. Damage initiation in each subcell is determined using quadratic 
definitions of damage strain eff e^, e % (a 3 -D extension of the strain-based Hashin criterion 17 ). 




where en are the normal strain components, 7 ij are the engineering shear strains, Q e , R e , and S e are the 
shear strain allowables, and the normal strain allowables X e , Y e , and are given by: 


Xj (Jii > 0 y _ ( F e T (J 22 >0 7 — \ ^ ^ 

If an<0 ’ 5 6 " \ Y? a 22 < 0 ’ 6 " \ Zf a 33 < 0 


( 5 ) 


where an are the normal stress components and the superscript L T\ or L C\ on the allowables indicates wether 
it is a compressive or tensile allowable. Damage initiates when any of the following inequalities are satisfied. 
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(6) 


III.B. Damage Evolution 

Once any of the inequalities in Eq. ( 6 ) is satisfied, damage begins to evolve in the matrix. It is assumed that 
the damage takes the form of an array of microcracks aligned perpendicular to the direction of the normal 
component of strain contributing to the damage strain (see Eq. ( 4 )) that has exceed one. These microcracks 
grow perpendicular to ay under mixed-mode (I, II, III) loading conditions, as shown in Figure 4 , where 
are the shear stress parallel to the microcrack faces. 
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Mode I Mode II 


Mode III 


Figure 4: Three loading modes that affect damage in the xi direction. 


As such, damage is introduced through six scalar damage variables Df and Dp, where i = 1,2,3 is 
determined from the damage strain(s) ef that exceeded one. It is assumed that the damage variables serve 
to degrade the engineering material properties. 


Ei = diE®; E 2 =d 2 E° 2 - E 3 = d 3 E° 

Vl2=di V 12 , V21=d 2 V2l\ 

^13 = ^1^13; ^31 = 4^31; 

^23 = d 2 V° 3 ; V 32 = d 3 U 32 


(7) 


where Ej are the degraded Young’s moduli, Ef are the virgin Young’s moduli, v l3 , and i/fj are the degraded 
and virgin Poisson’s ratios, respectively. The degradation factors, di are calculated from the scalar damage 
variables. 

di = < 1 > 0 ^ no sum on z, i = 1, 2, 3 (8) 

i 1 - b % D ? <*« < 0 

where and are parameters that allow the dependence of the moduli on the damage variables to be 
weighted differently in tension and compression. By default bJl = = 1. Note that, the Possion’s ratios are 

degraded in Eq. (7) in such a manner that the off diagonal terms of the compliance matrix remain constant 

Sij S°J, (/• / ./•). 

Microcrack evolution also affects the shear moduli. 


G 23 = (1 - b 42 D' 2 - b A3 Dl)G ° 2 3 

G 13 = (1 - b 51 Dl - b 53 D' 3 )G° 13 (9) 

G 12 = (1 - beiDl - b 62 D' 2 )G° 12 


where 


D* = 


Dj an > 0 
Df an < 0 


(10) 


Gij are the degraded shear moduli, and Gf are the undamaged shear moduli. The shear moduli are affected 
by a linear, weighted combination of two corresponding damage variables. The contributing damage variables 
are associated with the two normals of shear stress components related to the shear modulus of interest. 
The default weighting factors bij are equal to 0.5. 

Evolution of the damage variables Di is controlled by damage stress af versus damage strain ef consti- 
tutive laws, as shown in Figure 5. From this curve the damage increment can be determined. 


dfP 

dDi = (1 -Di- k[) 


(ii) 


where Di is the damage state at the previous increment, def is the damage strain increment, k[ = K/Ef 
and ki is the instantaneous tangent stiffness of the damage stress-damage strain curve. For full details on 
the derivation of Eq. (11), please refer to Ref. 10. In addition, Ref. 10 contains additional MMCDM theory 
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used to control the failure of the material due to meso/maco cracking using the fracture toughness of the 
material. The focus of this current work is limited to microdamage, so that portion of the theory has been 
purposely omitted. 



Figure 5: An incremental elastic damage stress versus damage strain curve. Note that X = V e , Y e , or Z e . 


IV. Schapery Theory for Microdamage Modeling at the Lamina Level 


ST is an anisotropic, thermodynamically-based, work potential theory for continuous fiber-reinforced 
laminae. 2,3,5 It is capable of predicting the non-linear response of a composite lamina with evolving mi- 
crodamage. When coupled with classical lamination theory, or the finite element method (FEM), ST can 
be used to predict the non-linear behavior of composite laminates, or structures. 6,7,15,18 ST was used to 
model several different laminates in Ref. 4; the predictions were extremely accurate up to the onset of more 
catastrophic mechanisms such as meso/macro-matrix cracking, fiber breakage and/or delamination. ST has 
been since expanded to account for these other mechanisms. 4,5,19 ” 22 However, the scope of this work is 
limited to microdamage evolution, so the original formulation of ST is used. 

ST assumes the total work potential Wt can be partitioned into an elastic, recoverable portion W and 
a dissipated portion S. 

Wt = W + S (12) 

Figure 6 shows a representative, uniaxial, non-linear stress-strain curve. The recoverable, elastic strain 
energy density W is shown as the area within the shaded triangle underneath a linear unloading line passing 
through the origin (ST assumes that there is no plastic deformation in the matrix, a reasonable assumptions 
for many PMCs 23,24 ). The dissipated potential S is the area under the stress-strain curve, above the 
unloading line. The internal state variable (ISV) S' is a measure of the amount of work potential used to 
facilitate structural changes in the material associated with microdamage. Microdamage is the nucleation 
and evolution of microcracks, shear bands, and fiber-matrix debonds in the matrix of the composite and is 
the primary mechanism for non-linearity in the stress strain response of a composite lamina up to the onset 
of more severe mechanisms such as transverse cracking, fiber breakage or delamination. 

Invoking the stationarity of the total work potential with respect to the ISV, and defining the strain 
energy density for a composite lamina yields the ST evolution equation for the microdamage ISV. 


1 ( 2 dE 22 , _2 dG 12 ^ 

2 \ e22 ~dS~ + ll2 ~dS~ ) 


= -1 


(13) 


The only elastic properties that remain in Eq. (13) are E 22 and G 12 because the other properties, used to 
define W , are unaffected by matrix microdamage. In addition to Eq. (13), microdamage healing must be 
precluded. 


^>0 


( 14 ) 
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Figure 6: Non-linear stress strain curve showing the recoverable, elastic potential W, and the dissipated 

potential S. 


For full details on the derivation of Eqs. (13)-(14), please see Refs. 6,7. 

The damaged transverse and shear moduli ( £22 and G 12 ) are related to their pristine values and the ISV 
through a set of microdamage functions e^, and gs . 


E22 


E° 22 e T s {S) a 22 > 0 . 
E° 22 e c s {S) cr 22 < 0 ’ 


G 12 = G° 12 g s (S ) 


(15) 


The microdamage curves are measured from three uniaxial coupon tests. In Ref. 4 numerous off-axis 
unidirectional and angle-ply laminates were tested, and all laminates produced the same gs(S) and es(S) 
{e$(S) or e^(S)) microdamage curves. The lay-up only affected the endpoint of the microdamage curve, 
which marks the transition to more severe damage and/or failure mechanisms. Eq. (15) can be used in 
Eq. (13) to solve for the ISV S for a given strain state. In turn, Eq. (15) can then be applied to calculate 
the degraded moduli. Finally, the damage moduli can be used to calculate the corresponding stress state. 
Note that, ST does not require an initiation criterion. For low levels of strain, S remains small, and due to 
the assumed microdamage functional form, does not affect the stiffnesses noticeably. However, as the strain 
levels increase, S is amplified, and the non-linearity arises naturally. 

The inputs to ST are the set of measured microdamage functions which relate the degradation in the 
engineering parameters to the ISV. The MMCDM (within GMC) was calibrated against the same experi- 
mental data used to obtain the ST inputs. In Section V.B the MMCDM is used to predict the stress-strain 
curves of the same laminates presented in Ref. 4 and compared to ST. The ultimate goal, presented in 
Section V.C, being to demonstrate the viability of GMC for ICME in a multiscale framework by varying 
the microstructure of the composite, at the lower scale, and observe its effects, at the higher scales, on the 
ST inputs and the performance of multi-angle laminates. The expectation is to achieve a similar level of 
accuracy to that of ST on the global scale while providing the ability to modify microstructural details at 
the constituent scale. Ideally, this scale linking could be extended further by tying the constitutive behavior 
of the matrix to a polymer processing model. Then, the inputs of GMC would be coming as outputs from 
other models and the constituent properties, composite microstructure and laminate lay-up could be tailored 
to meet application specific performance requirementes, thus achieving true ICME. 


V. Results 


V.A. MMCDM calibration 

Experimental data and ST results from Ref. 4 were used to calibrate and validate the GMC model. The 
elastic lamina properties measured for AS4/3502 are given in Table 1, where En is the modulus in the axial 
(fiber) direction, E 22 is the transverse modulus, z /12 is the Poisson’s ratio, and G \2 is the shear modulus. 
These properties were used to back out the elastic properties of the fiber and matrix, given in Table 2, used in 
the 7x7, doubly-periodic, square pakced RUC (shown in Figure 7) which represents continuous reinforcement 
in the yi -direction. The matrix (represented by green subcells in Figure 7) was assumed isotropic with a 
Young’s modulus E 171 and a Poisson’s ratio z/ m , and the fiber (shown in Figure 7 with blue subcells) was 
assumed transversely isotropic, where E{ x is the longitudinal modulus, E 22 is the transverse modulus, v { 2 
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is the axial Poisson’s ratio, and G 12 is the axial shear modulus. Ref. 4 reported a volume fraction of 58%. 
Note that even though the picture of the RUC in Figure 7 is 2-D, the subcell fields are full 3-D (they just 
do not vary in the yi -direction); so, for a transversely isotropic material, the transverse Poisson’s ratio 
is needed to fully characterize the stiffness tensor of the material. It was assumed that v{ 2 = ^ 23 - 

Table 1: Elastic properties for AS4/3502 lamina. 4 


Property 

Value 

En 

125.8 GPa 

E22 

9.31 GPa 

FL2 

0.329 

Gl 2 

5.10 GPa 


Table 2: Elastic properties for AS4 fiber and 3502 matrix constituents. 


Fiber Property 

Value 

Matrix Property 

Value 

E f 
^11 

214.3 GPa 

E 171 

3.53 GPa 

^22 

18.7 GPa 


0.35 

v 12 

0.32 



„/ 

0.32 



nf 

^12 

81.4 GPa 






P = l> -.ty 


Figure 7: RUC used in square-packed analyses containing 7x7 subcell architecture. Blue subcells are 

composed of the AS4 fiber material, and green subcells are composed of 3502 matrix material. 
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Data from three uniaxial coupon tests ( [zb 45 0 ]s, [15°], and [zb 30°]^ ) were used to characterize the 
three ST microdamage function inputs (e#), e^(5), and gs{S)) for AS4/3502. As mentioned before, many 
laminates were tested and produced the same microdamage curves. 4 These three laminates provided the 
most data before the onset of a more catastrophic mechanism. The experimental data from the same three 
coupon tests were used to calibrate the inputs to the GMC micromechanics model using the MMCDM for 
matrix damage. Following Ref. 10, k [ in Eq. (11) was taken to be an exponential function of the damage 
strain and two constants A and B. d 

k[ = Ae B> (16) 


where 



A T , e? > 0 .of B T , ef > 0 
A c , ef < 0 ’ { B c , ef < 0 


(17) 


All of the parameters calibrated for the MMCDM are presented in Table 3. 


Table 3: MMCDM parameters calibrated from [zb 45 0 ]s, [15], and [zb 30°]s experimental coupon data. 


Tensile Property 

Value 

Compressive Property 

Value 

XJ 

0.0125 

xf 

0.0287 

S, 

0.055 

- 

- 

A T 

1.0 

A c 

1.0 

b t 

1.0 

B c 

1.0 

bl 

0.75 

l f. 

ll 

0.1 

^4 ii ^5 it & Qi 

0.5 

~ 

- 


The resulting uniaxial stress strain curves for the [zb 45°]^, [15°], and [zb 30°] ^ laminates are presented in 
Figure 8. The MMCDM was able to correlate well for all three laminates. Some deviation can be observed 
at large strains in the [15°] and [zb 45°] ^ shown in Figures 8a and 8c, respectively. 

Figure 9 shows the microdamage curves measured from experiment 4 (that where also used as input for 
characterizing ST), and the microdamage curves generated from the output of the micromechanics model 
used to calculate the stress-strain curves in Figure 8. The degradation in the moduli (value of microdamage 
function) are plotted against a reduced ISV S r = S' 1 / 3 in Figure 9 (see Figure 6 for the definition of S). 
The use of the reduced ISV allows for the microdamage curves to be fit with polynomial functions. The 
micromechanics results correlate best with to e ^ for the [15°] laminate, shown in Figure 9a. 

The experimental microdamage curve for transverse (perpendicular to fiber) compression show an increase 
in moduli with increasing S r , see Figure 9b. There are pre-existing micro-cracks, -voids, and -fissures in 
the composite prior to being loaded. Some irrecoverable energy is required to close these flaws, and as 
they are closed, the transverse stiffness of the plies increases. After the flaws are fully closed, further 
compression acts to slightly reduce the stiffness; yet, the effect is minimal. Since a multiplicative damage 
variable is utilized in MMCDM an increase in stiffness would require a negative damage increment. However, 
inadmissability of healing, in MMCDM, is enforced on the damage increment. ST also uses a multiplicative 
damage variable, but this variable is equal to a function, and inadmissability of healing is enforced on the 
ISV argument of that function. Thus, ST can admit a positive change in stiffness with a positive change in 
the dissipated potential. To admit an increase in the transverse stiffness with MMCDM, the “no healing” 
restriction imposed on the damage increment would need to be relaxed. This is physical incorrect though 
because damage is still increasing (i.e. energy is dissipated to advance structural changes in the material); 
under transverse compressive loading, the dissipated potential serves to advance changes that manifest as an 
increase in the observed stiffness. The discrepancy, between MMCDM and the experiment, in the transverse 
behavior of the plies had little effect however on the stress-strain response of the [zb 30°] ^ laminate, as seen 
in Figure 8b. 

The shear microdamage curve obtained from MMCDM deviates from the experimental curve for low 
values of S r (Figure 9c). This is because MMCDM require an initiation criterion, and G 12 /G 12 will remain 
one until initiation criterion is satisfied. In the experiment, initial, low levels of dissipated potential act to 
degrade the shear modulus; however the degradation is small and does not significantly impact the stress- 
strain response of the [zb 45°] ^ - As S r increases the deterioration of G 12 is more apparent. ST can capture 
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Figure 8: Uniaxial stress-strain curves used to calibrate MMCDM and characterize ST. 


11 of 21 


American Institute of Aeronautics and Astronautics 




S r (MPa 1 / 3 ) 


(b) Transverse Compression. 



Figure 9: Matrix microdamage curves from experiment and GMC using MMCDM. 
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this behavior directly because it does not require a failure criterion. Although there is a noticeable difference 
in the shear microdamage curve from MMCDM, compared to experiment, it is not evident in the stress-strain 
curves shown in Figure 8c. However, this indicates that the functional form of the damage evolution in the 
MMCDM, Eq. (11) may not be correct. If Eq. (11) utilized based on work potential, rather than strain, 
perhaps the predictions of the curves in Figure 9 using MMCDM would be improved 

V.B. MMCDM Predictions 

The calibrated properties, in Tables 2 and 3, were used in the MMCDM within GMC to predict the response 
of several different off-axis unidirectional and angle-ply laminates. The resulting stress versus strain curves 
are compared to experimental results from Ref. 4 in Figure 10. The focus of this work was to evaluate 
the capability of the MMCDM to predict the onset and evolution of matrix microdamage with the intent of 
subsequently simulating changes in the composite processing and manufacturing to improve the performance 
of the material, with respect to microdamage. Many of the experiments used to produce results presented 
in Figure 10 exhibited damage mechanisms beyond matrix microdamage, including transverse cracking and 
delamination. ST was validated as an accurate tool for predicting microdamage in PMC laminates, 4 but de- 
viated from the experimental results upon the onset of other damage mechanisms. Thus, the ST predictions, 
taken from Ref., 4 are also included in Figure 10 to give a better indication of the ability of the MMCDM 
to predict microdamage. However, MMCDM is not restricted to modeling microdamage alone, and can be 
used to model the combining effects of all damage mechanisms. Whereas, ST would require a reformulation 
to capture the effects of other mechanisms. 

The GMC model utilizing MMCDM is able to predict the non-linear response of simple laminates ex- 
hibiting matrix microdamage reasonably. Figures lOa-lOc show some deviation from the ST predictions and 
experimental results at high levels of strain, especially in the [45°] and [60°] laminates. For the angle-ply 
laminates (Figures lOd - lOg), the GMC/MMCDM model was able to predict the non-linear behavior of 
the [zb 15°]^, [A SO 0 ]^, and [d= 60°]s laminates well, but there is significant error in the [=b 40°]^ laminate 
for strains greater than 2.5%. The [zb 40°] ^ laminate was loaded to four times the strain level of the other 
angle-ply laminates. As such, the disparity observed at the end of the [zb 45°]s, at approximately 3% strain, 
in Figure 8c continues to amplify with increasing strain in the [zb 40°] ^ example. This suggests that a differ- 
ent functional form may need to be introduced into the MMCDM to more adequately capture the damage 
behavior. 

Figure 11 shows the lamina microdamage curves predicted from the simulation of the same laminates, 
used to generate the stress-strain curves in Figure 10, with GMC utilizing the MMCDM. These microdamage 
curves are input for the ST macroscopic damage model, and as inputs, do not change based upon lay-up. Ref. 
4 reported that these microdamage curves were also insensitive to changes in the lay-up experimentally, with 
the exception of the endpoint of the microdamage curves which marks the initiation of some other damage 
or failure mechanism. Thus, only a single curve is presented for each microdamage function, representing 
the experimental data. 

The transverse tensile microdamage functions in Figure 11a predicted using GMC with MMCDM exhibit 
a fair degree of scatter. However, this did not directly correlate to poor accuracy when predicting the 
stress-strain curves, shown in Figure 10. The transverse compression microdamage curves are more or less 
coincident, but there is a an appreciable difference between the predicted curves and the experimental data 
for reasons explained previously. The predicted shear microdamage curves also exhibit accordance with each 
other in Figure lib. 

The laminates were loaded, in strain control, to the maximum strain reported in Ref. 4. At this point 
the laminates either failed or exhibited some other catastrophic damage mechanism. The endpoints of the 
predicted microdamage curves are marked with red symbols in Figure 11 and demonstrate that the transition 
to failure, or other damage mechanisms, is dependent on the lay-up and does not occur at a fixed value of 
S r . This means that the driving force for microdamage is independent of the lay-up, but the driving force 
for other mechanisms is not. 

Figure 12 shows the stress-strain curve predictions for three multi-angle laminates: [±45°/ zb45°/90°/ zb 
45 °] s , [±45°/90§/ ± 45°/90£/ zb 45°] 5 , and [±45°/60£/ zb 45°/60£/ ± 45°] 5 . Reasonable predictions were 
obtained with MMCDM for all three laminates, and these predictions match those for ST. Ref. 4 reported 
the initiation of transverse cracks at axial strains of 0.004-0.005; however, the ST results do not deviate 
noticeably from the experiment until larger strains are achieved. This indicates that, initially, the effects of 
a small number of transverse cracks on the laminate stiffness is minimal but becomes more prominent as the 
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(a) [30°]. 


(b) [45°]. 




(c) [60°]. 


(d) [±15°] s . 




(e) [±40°] s . (f) [±50°] s . 



(g) [±60°] s . 


Figure 10: Stress-strain curves predictions from GMC with MMCDM compared to experimental results 4 

for a variety of off-axis unidirectional and angle-ply laminates. 
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(a) Transverse. 



(b) Shear. 

Figure 11: Matrix microdamage curves from GMC using MMCDM for a variety of laminates. 
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crack multiply. 




(a) [±45°/ ± 45 o /907 ± 45°] s . (b) [d=45°/90§/ ± 45°/90§/ db 45°] s . 



(c) [±45760°/ ±457607±45 o ]s. 


Figure 12: Stress-strain curves predictions from GMC with MMCDM compared to experimental results 4 

for a three multi-angle laminates. 


Figure 11 clearly indicates that the current damage model used to represent the matrix damage within 
GMC does not produce a commensurate amount of change in global stiffness with respect to the ISV S 
(measure of dissipated work), nor are the predicted microdamage curves insensitive to laminate lay-up, 
as observed in the experiments. However, ST is able to predict the stress-strain behavior of composites 
experiencing matrix microdamage because the ISV chosen yields the correct driving thermodynamic driving 
force / for microdamage evolution. The thermodynamic driving force is given by the change in the elastic 
strain energy density with resect to incremental changes in the ISV. 3 


dW 

~dS 


(18) 


If S was chosen inappropriately, then the thermodynamic driving force in the model would be incorrect. 
In addition, the functional forms of the degradation of £22 and G \2 were also chosen correctly (measured 
directly from experiment), and moreover, are insensitive to the lay-up of the laminate. 

The MMCDM is not a work potential theory, and accordingly does not utilize a thermodynamic force 
to control the evolution of microdamage. Yet, the parameters can be tuned to mimic the characterization 
experiments for ST and predict the response of other laminates within a fair degree of accuracy. Unfortu- 
nately, this makes it difficult to assess wether it is the constitutive model, the micromechanics, or both that 
are responsible for the observed discrepancies. To reduce potential error induced by the constitutive model 
a 3-D work potential theory for initially isotropic materials 15 will be used in the future to model the matrix. 
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Confidence in the constitutive law could be achieved if the predicted microdamage curves were insensitive to 
the laminate lay-up, as observed in the experiments. It should also be mentioned that the flexibility of the 
MMCDM allows multiple mechanisms to be captured (although not demonstrated in this work), whereas, ST 
would require the addition of another ISV 4,5,19 significantly increasing the complexity and computational 
cost of the theory. Moreover, carefully planned experiments must be devised to isolate and/amplify the 
effects of the mechanism of interest, relative to other mechanisms, and it must be determined that the choice 
of ISV for the new mechanism is correct and independent of lay-up and loading. 

V.C. Using GMC and MMCDM to Predict Processing Effects 

Although using the MMCDM within GMC was not able to identically reproduce the macroscopic, ST results, 
it can still be employed to interrogate the effects of microstructures which is linked with processing, such 
as fiber volume fraction and fiber packing, on the microdamage evolution in the matrix phase. Thus, 
demonstrating the potential usefullness for GMC within ICME applications. This type of study is not 
possible with ST, or other macroscopic theories. 

The fiber volume fraction of the [15°], [TSO 0 ]^, [±45 0 ]£ were varied from 10% to 78.5% assuming the 
square packed architecture shown in Figure 7 to determine the effects of fiber volume fraction on the ST 
microdamage functions. 

Figure 13 shows the normalized degradation in the appropriate moduli versus fiber volume fraction V/, for 
a fixed damage state (S r = 0.5), obtained from simulations of the three coupons. The transverse compression 
response is completely insensitive to changes in Vf. Under transverse tension the degradation of E22 increased 
by 5% as Vf increased. However, the damage began to lessen as Vf goes above 60%. The shear response 
was extremely sensitive to the fiber volume fraction and the damaged shear modulus steadily declined with 
an increase in Vf. After Vf = 50% the rate of degradation with respect to increasing Vf became very rapid. 
According to Figure 13, a lower fiber volume fraction will lessen the severity of any resulting damage on 
shear modulus of a ply. Although, there will be a penalty in stiffness, the fixed damage state (S r = 0.5) may 
be achieved at lower loads; therefore, there must be a balance. 



Figure 13: 


Effect of fiber volume fraction on modulus degradation for a fixed damage state (S r 


0.5). 


Since Vf affects the stiffness of the laminates, comparing stress-strain curves does not accentuate the 
effects of fiber volume fraction on the microdamage evolution in the composite. Therefore, Figure 14 shows 
the fraction of the total work that is dissipated due to microdamage, for a fixed total work potential Wt = 2.5 
MPa, as a function of fiber volume fraction, for the three characterization laminates and the three multi- 
angle laminates. Most of laminates exhibited an increase in dissipated potential with fiber volume fraction. 
The amount of potential that was dissipated is more sensitive to changes in the Vf for the [±30°]^ and 
[±45°] s' angle-ply laminates. The second multi-angle laminate, which had 40% 90° plies, actually showed 
slight decrease in dissipated potential up to 40% Vf followed by an increase in S/Wt • Figure 14 shows that 
a lay-up and volume fraction can be chosen to minimize the fraction of the total, applied work potential that 


17 of 21 


American Institute of Aeronautics and Astronautics 


is used to facilitate microdamage evolution, but again, attention must be paid to stiffness and load carrying 
requirements of the structure. 



Figure 14: Fraction of total work potential that is dissipated versus fiber volume fraction for a fixed total 
work potential Wt = 2.5 MPa. 


A preliminary investigation on the effect of packing on the performance of the material was conducted 
by changing the fiber architecture from a square packing (Figure 7) to a hexagonal packing, shown in Figure 
15. The volume fraction was retained at 58%, as well as the fiber geometry. 


Vf=5..fl00GQQe-G1 



*2 


Figure 15: Hexagonally-packed fiber architecture. 


Figure 16 shows the predicted microdamage curves for the two different microstructures using the charac- 
terization laminates ( [15°] , [TSO 0 ]#, [±45°]^). The transverse tension and shear microdamage functions were 
only moderately sensitive to changes in the packing, but the transverse compression microdamage function 
was significantly affected by the change. However, it cannot be stated conclusively that e$ is sensitive to the 
fiber packing because a good correlation for e$ was not obtained for the baseline, square-packed architecture. 

The stress-strain curves for the three multi-angle laminates predicted using MMCDM within GMC uti- 
lizing square and hexagonal packing architectures are displayed in Figure 17. The packing appears to have 
a minimal effect on the results for the multi-angle laminates. However, based on the results presented in 
Figure 16, the difference is expected to be more pronounced if the laminates are loaded in compression. 
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Figure 16: Matrix microdamage curves from GMC using MMCDM for two different fiber packing architec- 
tures. 



€-XX 

Figure IT: Stress-strain response of three multi- angle laminates for square-packed and hexagonally-packed 
architectures. 
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It is difficult to make any definitive statements on the effects of fiber packing on the microdamage in 
the composite with only two cases. To make any solid conclusions the statistical repeating volume element 
(RVE) containing randomly packed would first need to be identified and used as a baseline for comparison 
with other perturbations of this system (the extreme cases being ordered microstructures such as square and 
hexagonally packed), as having the correct statistical RVE is important with respect to not only accurately 
concluding about packing effects but also whether or not one has the correct damage model. Moreover, 
the need to have accurate constitutive models are imperative, more so with micromechanics than at the 
macroscale since the in-situ fields are highly multiaxial and nonproportional. 

VI. Conclusion 

The MMCDM was implemented within the MAC / GMC software suite and used to predict the response 
of numerous off-axis unidirectional, angle-ply, and multi-angle laminates exhibiting matrix microdamage. 
The matrix damage parameters used in the model were calibrated to match the response of three different 
coupon tests. These same three coupons were used in Ref. 4 to characterize the microdamage functions 
used in macroscopic ST damage model. The stress- strain curves predicted using MMCDM were compared 
to the experimental data and ST results. In addition, the microdamage functions predicted by MMCDM 
were compared to experiment. 

Figures 8-12 demonstrate that the MMCDM, within a micromechanics theory, can correlate very well 
to data used to characterize the ST macromechanics model. However, the MMCDM did perform poorly in 
predicting the response of some laminates and exhibited a marginal amount of error when reproducing the 
matrix microdamage ST input functions. The discrepancy between ST and MMCDM was less apparent in 
multi- angle laminates. Ideally, an work potential-based damage theory, similar to ST, should be introduced 
into MAC/GMC to model the behavior of the matrix constituent. 

ST offered accurate predictions for the response of laminates experiencing matrix microdamage evolution 
making it an effective constitutive model for use in FEM modeling of structures. However, the influence of 
changes in the material processing or microstructure cannot be determined with ST. With micromechanics 
(GMC), and MMCDM, parameters such as volume fraction and fiber packing can be traded and their effects 
on the damage evolution of the system can be observed. Furthermore, the micromechanics can be utilized 
to predict the ST microdamage functions, accounting for processing changes. 

It was seen that changes in the volume fraction significantly affected the shear microdamage function. In 
addition, the fraction of total work potential that was dissipated due to microdamage (for a given applied 
total work potential) was sensitive to changes in the volume fraction for particular laminates, but insensitive 
to others. Finally, it seemed that changing the fiber packing architecture from square to hexagonal has 
a marked impact on the resulting transverse compression microdamage function. Further studies can be 
conducted to determine the effects of random architectures on the response of the laminates, but this will 
require a thorough statistical investigation to identify the statistical RVE. It is anticipated that the damage 
curves produced using a statistical RVE, within a laminate, will not be sensitive to the lay-up of the laminates 
because the statistical RVE is energetically objective when damage is present. 25,26 

Through these simulations, it has been demonstrated that GMC can be used as a viable tool for conduct- 
ing ICME on composite materials while achieving comparable macrolevel results to an accurate macrolevel 
theory ST which lacks the ability to design the material. However, these are preliminary steps towards ICME, 
but to achieve true ICME the laminate stiffness, load-carrying capability, and susceptibility to microdamage 
must be tailored by varying the fiber volume fraction, microstructure, laminate lay-up, and material prop- 
erties to meet some structural performance. Additionally, the behavior of the constituents should be tied to 
more fundamental materials science models to improve the robustness of the method and its applicability 
to ICME. Moreover, the properties that are determined suitable to meet the required performance must be 
tied to material manufacturing and processing simulations or data to determine if they are achievable. 
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